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The insulating ternary oxide YbFe2 04 displays an unusual frustration-driven incommensurate 
charge-ordering (CO) transition, linked to possible ferroelectricity. Based on high-resolution syn- 
chrotron data, we report a detailed structural model showing that the CO phase is an incommen- 
surate charge- density wave and cannot be ferroelectric, since the electrical dipole moments are also 
incommensurately modulated. The change between continuous and "spotty" helices of scattering 
at the CO transition is attributed to three-dimensional fluctuations of the direction of the ordering 
wavevector. 



PACS numbers: 77.80.-e, 71.45.Lr, 61.05.cp, 75.85.+t 
I. INTRODUCTION 



The family of systems with general formula RFe204 
{R = Dy-Lu, and Y) has attracted significant interest, 
due to the possibility of high-temperature multiferroic 
behaviour arising from charge orderin^^^n^. Just above 
room temperature these systems undergo a metal-to- 
insulator (MI) transition, associated with sharp anoma- 
lies in the specific heat^, and on further cooling magnetic 
ordering sets in at around 250 k!^31III^ xhe MI tran- 
sition has been interpreted as a peculiar form of charge 
ordering (CO), where the mixed- valence character of iron 
(the average Fe valence is 2.5+) competes with the ge- 
ometrical frustration of the lattice. The crystal struc- 
ture of RFe204^^ consists of an alternating stacking of 
RO2 layers and Fe2 02 bilayers (with three such bilayers 
per conventional unit cell). The Fe ions form a trian- 
gular pattern within the layers, so the system is geo- 
metrically frustrated with respect to both the Fe spins 
and valenc^^ and it is impossible to order an equal 
amount of Fe^+_and Fe^+ on the same layer. Yamada 
and coworker^I^'I^ proposed a possible CO structure for 
LuFe204, which they predicted to be an incommensurate 
charge-density wave (CDW) and an incommensurate an- 
tiferroelectric. Although their work is mainly based on a 
theoretical model (see also below), they were able to ex- 
plain the incommensurate positions of the CO satellites 
which, in reciprocal space sections, 'meander' around the 
lines of commensurability. Ikeda and coworkers later 
reported that the low-temperature insulating phase of 
LuFe204 is ferroelectric^, and proposed a much simpler 
commensurate ferroelectric CO model in which two in- 
equivalent lattices would be formed within the bilayer: 
Fe^++2Fe^^+^^+ on one layer and Fe^++2Fe^^"^^+ on the 
adjacent one (Fig. [l^), creating net dipole moments for 
each bilayer, which would stack ferroelectrically. More 
recently, Angst et a/.^ presented an alternative model, 
based on LuFe204 diffraction data and a commensu- 



rate approximation of the charge-ordered superstructure, 
and proposed that the bilayer stacking is antiferroelectic 
rather than ferroelectric. 

In this paper, we present a set of three-dimensional 
reciprocal space data on YbFe204. We focus on the 
R = Yb system since we found it to have more sta- 
ble oxygen content (t o whi ch the properties of the sys- 
tem are very sensitivd^^'^ and to yield better-quality 
single crystals than R = Lu. These data clearly show 
that the meandering lines of incommensurate satellites 
are actually helices, as predicted theoreticall}'^^ We 
also present a real-space model of the fully incommen- 
surate CO structure, which, in spite of its simplicity, is 
in very good agreement with the experimental data. We 
show that, below T* ^ 320 K, YbFe204 forms a true 
incommensurate CDW (as predicted by Yamada) asso- 
ciated with an incommensurate rotating strain pattern. 
Although local electric dipole moments are present, their 
directions reverse periodically, implying that each bilayer 
has a zero net moment. Therefore, the previously ob- 
served electrical and dielectric anomalies are most likely 
associated with CDW physics rather than with ferroelec- 
tricity. Above T* the individual superlattice peaks melt 
into continuous helices of scattering, which nonetheless 
remain very coherent. This may indicate that low- lying 
incommensurate states with different wavevectors are be- 
coming progressively populated — a situation that is not 
observed in conventional CO transitions. 



II. EXPERIMENTAL METHODS 

YbFe204 single crystals were grown by the opti- 
cal floating-zone technique. Stoichiometric mixtures of 
Yb203 and Fe203 were ground together and calcined 
in C0:C02 (4:5) mixed gas flow at 1160°C for 48 h 
with intermediate grinding. The feed and seed rods (10 
mm diameter) were finally sintered at 1175°C for 36 h 
in C0:C02 (1^2) mixed gas flow. The crystal growth 
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was carried out in C0:C02 mixed (1:4) atmospheric flow 
using a IR-heated image furnace (Crystal System Inc.). 
The feed and seed rods were rotated in opposite direc- 
tions at 25 rpm during crystal growth at a rate of ~2 
mm/h. Preliminary X-ray diffraction data were collected 
on an Agilent Technologies Supernova diffractometer, us- 
ing a Mo source (wavelength 0.71073 A) and cryostream 
to cool to 150 K. Synchrotron X-ray diffraction exper- 
iments were carried out on beamline 119 at the Dia- 
mond Light Source, UK. The beamline employs a double 
bounce monochromator and focussing mirrors to bring 
the beam down to approximately 170 /im x 85 /im at 
the sample position. A wavelength of 0.6889 A was used 
throughout and data were collected using a Rigaku Sat- 
urn 724 CCD detector at temperatures between 150 and 
360 K. Reciprocal space sections were generated from 
the raw CCD images using the programme CrysAlis^^°. 
Stacks of such sections were converted to a volume data 
set and analysed using FIJl^^ to enable visualisation of 
the scattering intensities in three-dimensional reciprocal 
space. 



III. RESULTS 

Figs. [l]3-[l]i show reciprocal space sections of X-ray 
diffraction data collected from both the laboratory source 
and synchrotron. The (/i, /c, 12) section (data at 150 K 
are shown in Fig. [TJd; the same section at 360 K looks 
very similar) clearly shows the appearance of satellite 
spots, which can be indexed approximately on a x ^/3 
supercell. A perpendicular section through these satel- 
lites (Fig. ^) shows that the scattering is densely dis- 
tributed around (m + 0.666, —2.33, /) lines, with m an 
integer. The Bragg peak maxima are not at commensu- 
rate positions, but oscillate in h and k around these lines 
and occur at incommensurate values of separated by 
- 1/3 c*. Above T* (Fig. [l]i), the Bragg peaks 'melt' 
into a continuous meandering line (this intensity is re- 
ported to vanish above ~ 500 Kf^ where CO is no longer 
present), but remain very coherent perpendicular to the 
line itself. The data in these sections are very reminiscent 
of previously published LuFe204 data (for example ref. 
H]), although in LuFe204 the spacing along the / direc- 
tion is 1/2 c*. However, our high-resolution synchrotron 
measurements were suitable for full three-dimensional 
reconstructions (rather than simple sections), affording 
much richer views of the experimental data. The three- 
dimensional iso-intensity contours of scattering for one 
rod (Fig. [l^) and for a larger region of reciprocal space 
(Fig. 2| evidence that the meandering lines in Figs. ^ 
and[Ti are nothing but sections through helices of scat- 
tering in reciprocal space. The radius of the helix, p, 
which is of the order 0.025a* at 150 K (and approxi- 
mately twice this value at 360 K), indicates that the 
^/3 X \/3 supercell is only a commensurate approxima- 
tion to a complex incommensurate CO superstructure. 
Figs. [2^ and show a three-dimensional view of the 



'spotty helices' below T* and of the continuous helices 
above T*. 

These deviations from the commensurate structure 
arise because the commensurate arrangement of charges 
proposed by Ikeda et al^ is energetically unstable due 
to the frustrated topology of the rhombohedrally-stacked 
bilayers. Rastelli and TassP^ flrst explored the effects 
of this frustration in the context of magnetic ordering 
in solid oxygen, which displays a topologically analogous 
stacking of magnetic triangular layers, and found that the 
'conventional' 120° magnetic ordering, with propagation 
vector (1/3, 1/3, 0), is unstable. Instead, if interactions 
between the planes are taken into account, the system 
displays a line of degenerate ground states, correspond- 
ing to reciprocal space q vectors lying on the contour of 
a helix that surrounds the ( 1/3, 1/3, 0) point. In his 
previously cited analysi^^-^^, Yamada et al. exploited 
the topological similarity between the i?Fe204 family and 



o • O • • \ O 
O* •* •* o 

a Te(3-sH 





(d) 360 K 



(/z,-2.33,/) 



h [r.l.u.] 



(e) 




(0.33,-2.66,/) 



FIG. 1: (Color online) (a) The simplest ordering model of 
Fe^^ and Fe^^ on adjacent FeO layers in i?Fe2 04^- Lower 
panel: oxygen displacement pattern in a single FeO plane. 
The arrows can also be used to illustrate the conventional 
120° magnetic structure of a triangular Heisenberg magnet. 
The CO and crystallographic unit cells are indicated, (b) 
& (c) Reciprocal space sections from lab data collected on 
YbFe204 at 150 K, respectively the (/i, k, 12) and (h, -2.33, /) 
layers. Bragg peaks in (b) are subject to the selection rule 
—h-\-k-\-l = 3n and are surrounded by hexagons of satellites, 
(d) Synchrotron data showing the (/i, —2.33, /) layer at 360 
K: the streaks observed at low T (b) become more diffuse for 
T > T* . The satellites shown in (c) and (d) are marked with 
a box in (b). (e) Three-dimensional iso-intensity plot of the 
(0.33, —2.66,/) satellite at 150 K (synchrotron data). 
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the magnetic systems considered by Rastelli and Tassi, 
and modelled the 'wavy' pattern of diffuse scattering and 
Bragg peaks observed in LuFe204 in terms of recipro- 
cal space 'helices'. More recently, Harris and Yildirim^^ 
constructed a complete mean-field theory of the CO in 
i?Fe2 04, and similarly to the magnetic case they con- 
cluded that the mean-field ground state (calculated us- 
ing interactions up to the third neighbour) results in a 
degenerate helix of q- vectors. A fourth-neighbour term 
would stabilise a single ground state, which should dis- 
play Bragg peaks at h — l/3-h^i,A^ = 1/3+^2 and / either 
integer or half-integer depending on the sign of the inter- 
action. The mean-field models in refs. 18a and ,21, are best 
compared to our high-temperature data, since their fully 
ordered phases (with four interactions) have Bragg peaks 
in different places along the helix. 

The topology (handedness, phase, and pitch) of each 
satellite helix is fixed entirely by the (chiral) holohedral 
axes that are present along c* in the rhombohedral re- 
ciprocal lattice at (1/3,1/3) type positions. Thus, for a 



reciprocal lattice of points satisfying —h-\-k-\-l = 3n 
where n G Z, the chirality of each helix is given by 
h -\- k = (3n ± l)/3 where + (— ) corresponds to left- 
(right-) handed helices (Fig. |2]i), and the relative phase, 
Acp^ of adjacent helices must be 27r/3 (Fig. [2|c)). The 
pitch of the reciprocal space screw axes additionally fixes 
the pitch of the satellite helices to 3c*, whilst the radius 
of the helices is an adjustable parameter in the theory, 
and corresponds to the ratio of the nearest- and next- 
nearest-neighbour interactions. The 150 K data (Fig. ^) 
indicate that at low T unique modes have been stabilised 
with peaks appearing along the helices with characteris- 
tic separation of ~ 1/3 c*. Although this is at odds with 
the theory, the mean-field model could be modified, in- 
troducing longer-range interactions, to yield the correct 
q-vector. The intensities of the peaks also follow the 'en- 
velope' of the high T data, indicating that the structural 
distortions involved are very similar in both phases. 



IV. CONSTRUCTION OF THE 
INCOMMENSURATE MODEL 
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FIG. 2: (Color online) Three-dimensional reciprocal space iso- 
intensity plots of the satellite helices, (a) and (b): Helices 
below and above T* . At low T intensity is concentrated onto 
peaks separated by 1/3 c*. (c) Three adjacent helices (T = 
360 K). The pitch of each helix is 3c*, but the phase varies 
according to the position in the (a* , 6*) plane, (d) Two helices 
at T = 360 K, for 24 < / < 10 r.l.u., illustrating the selection 
rule for the handedness (denoted RH/LH for right/left), (e) 
Schematic reciprocal space reconstruction of several satellites 
at 150 K. The reciprocal lattice is shown below (white circles) 
with the locations of the satellites (red/dark circles) projected 
downwards. 



While our analysis supports Yamada's original pro- 
posal, conclusive evidence can only come from a struc- 
tural model that reproduces not only the Bragg peak 
positions, but also the observed intensities. Although 
the present data cannot be reliably integrated, prevent- 
ing a full crystallographic refinement, we present an ex- 
tremely simple, parameter- free model of the incommen- 
surate CDW that yields very good agreement with the 
experimental intensity profiles. 

Initially we construct a commensurate displacement 
(strain) pattern of the oxygen ligands that is consis- 



O 
i 

c 

Fe2+ 



Yb 







o 



?4 







o 

0^ • .^o 



o 



. Vertical 
displacement 



In-plane 
displacement 



Vertical 
displacement 



a 



FIG. 3: (Color online) Schematic picture of the oxygen dis- 
placements within the YbFe204 structure. Apical oxygens 
(within the Yb02 layers) are displaced vertically according 
to a sinusoidal modulation (dashed line), whilst the oxygens 
in the Fe202 layers are displaced in plane (part of the dis- 
placement pattern is shown by the arrows). As shown, the 
modulations follow the commensurate v^x structure. The 
same modes are propagated with incommensurate q to gen- 
erate the final structure. 
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tent with the CO structure proposed by Ikeda et al^. 
Such patterns, where metal oxygen distances are smaller 
(larger) for higher (lower) oxidation states respectively, 
are ubiquitous in CO metal oxides. Ikeda's CO structure 
can be reproduced by a single commensurate wavevector 
q = (1/3, 1/3,0) and two displacements for the in-plane 
and apical oxygens (along the fixed (210) and [001] direc- 
tions, respectively: see Fig. [3|. The oxygen displacement 
pattern in the Fe2 02 layers is completely defined by this 
rule and by the requirement that the two layers be in 
phase opposition. This yields six kinds of distorted pyra- 
mids, shown in Fig. [4^ with the corresponding valence 
states. A further constraint can be imposed by assuming 
that the in-plane displacement of each oxygen atom is 
proportional to the value at its position of a 'distortion 
field' S obeying the relation 



VS 



-5p 



(1) 



where Sp is the charge (valence) density deviation from 
the average (this yields essentially equal displacements on 
the two oxygen sites and a one-parameter model). Eqn. 
([T]) is reminiscent of Maxwell's equation for the electric 
field, but here the displacements are of steric rather than 
electrostatic origin. The in-plane displacement pattern 
is the polar analogue of the 'conventional' 120° magnetic 
structure: therefore, at a deeper level, eqn. ([T]) defines 
a relation between scalar and vector ordering on a frus- 
trated lattice - a concept that underpins, for example, 
magnetic monopoles in spin ic^^. A minimal micro- 
scopic model of the incommensurate phase is then ob- 
tained simply by employing the correct incommensurate 
wavevector, producing, for the in-plane oxygens, the po- 
lar analogue of the incommensurate magnetic structure 




FIG. 4: (Color online) Real space diagrams of YbFe204- (a) 
The oxygen distortions involved in the six possible pyramids 
in the x structure, (b) Example of an incommensurate 
oxygen displacement pattern in the FeO planes, made up from 
a continuous variation of the configurations in (a). The irons 
are shaded according to their effective valence determined by 
eqn. ([T]), and the direction of q = (1/3+^, 1/3+^) is indicated 
by a bold arrow. 



proposed in ref. [23l This model yields a continuous vari- 
ation of the Fe valence on the different sites (with a cor- 
responding variation between the configurations in Fig. 
|4^). This is shown pictorially in Fig. ^p. The magnitude 
of the in-plane and out-of-plane oxygen displacements is 
assumed to be equal for simplicity. 



V. COMPARISON WITH SCATTERING DATA 

The experimental data at 150 K show that the scatter- 
ing intensity is located at positions along the rod sepa- 
rated by an incommensurate value of Qz ~ 1/3. Such an 
intensity distribution cannot be reproduced with a sin- 
gle wave vector, but requires two symmetry- inequivalent 
ordering wave vectors (propagation vectors) and their 
symmetry-equivalent domains as shown in Fig. |5j The 
first, qi = (1/3 + (^1, 1/3 + 62^ 1/3), marked with (*) in 
Fig. |5^, gives rise to 12 equivalent propagation vectors of 
which six (solid arrows in the figure) contribute to each 
satellite rod {Si > 62). Considering a structure with only 
this propagation vector would lead to peaks at / = n±l/3 
for integer n, positioned around the helix. In addition, to 
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FIG. 5: (Color online) (a) Reciprocal space sketch indicat- 
ing contributions from the two propagation vectors, qi and 
q2, to the helical scattering pattern. The black circles rep- 
resent Bragg peaks and the grey shows the location of the 
satellite helix, (b) & (c) Calculated intensities along the 
(1/3,-2/3,/) rod from the qi and q2 structures respectively, 
(d) Overall pattern expected given equal contributions from 
both structures, which is compared against the data in Fig. 
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FIG. 6: (Color online) Reciprocal space sections and simulations at T = 150 K. Data (left-hand panels) and calculations 
(right-hand panels) are shown for six reciprocal space cuts. The calculated intensities are derived from the oxygen displacement 
model described in the text. The intensities (all of which are scaled to those in (a) and (b)) are indicated by the color bars. 



obtain scattering at / = n we take a second propagation 
vector, q2 = (1/3 + e, 1/3 + e, 0) which results in a fur- 
ther three equivalent propagation vectors (dashed arrows 
in the figure). It seems highly plausible that if the two 
modes are closely separated in energy then the system 
may order in such a way that both modes are populated 
at 150 K. 

Scattering intensities due to the two qs were calcu- 
lated separately as follows. The radius of the helix, p ^ 
0.025a*, was used in order to establish commensurate ap- 
proximations to the incommensurate wave vectors: these 
were qi = (0.336,0.304,1/3) and qs = (0.320,0.320,0). 
For the qi structure, a supercell corresponding to this ap- 
proximate structure was generated and the oxygen atoms 
displaced according to the model described above. The 
oxygen displacements are a measure of the CO, which 
could be estimated, for example, using the bond valence 
sum method^, provided that both main peaks and satel- 
lites are measured. Here, only satellites were measured, 
and the displacements were nominally taken to be 0.1 A. 

The structure factors were calculated following the 
standard expression F{Q) = ^ fj exp((5a^r j + Qy^l + 
Qz^i) where the summation runs over all j atoms in the 



supercell and the fj are the atomic scattering factors 
(which depend on |Q|). The positions of the Yb ions are 
shifted slightly so as to be at the center of mass of the 
surrounding oxygens (imposing the periodic boundary 
conditions associated with the supercell). The resulting 
intensities (calculated separately for each domain) were 
then multiplied by an array of two-dimensional Gaussian 
peaks at the positions given by the propagation vectors 
to yield the predicted scattering pattern (shown for the 
(1/3, -2/3, 1) satellite helix in Fig. [5]d). The width of this 
Gaussian reflects the coherence lengths of the wave vec- 
tor domains (~ 150 A at 150 K). Intensities from the q2 
structure were calculated in an identical way (Fig. [sj)). 
In order to best model the data we assume that an equal 
population of qi and q2 modes exist in the crystal and 
therefore simply average over the two intensities to yield 
the final pattern (Fig. |5]i). Given that there is no ad- 
justable parameter other than a single scale factor, the 
observed scattered intensities along the helix are mod- 
elled with surprisingly good agreement over the whole 
reciprocal space (see Fig. [6|, although some of the more 
subtle features are slightly different, as is to be expected 
from such a simple model. 
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VI. DISCUSSION AND CONCLUSION 

Our analysis of high-resolution synchrotron X-ray scat- 
tering data collected on YbFe204 both above and below 
T* provides new insight into the physics of this class of 
materials. Below T*, we can model the helical envelope 
and intensity distribution of the satellite peaks in recipro- 
cal space using an incommensurate displacement model. 
Since the stoichiometry would favour a simple y/3 x ^/S 
ordering, our observations imply the presence of an in- 
trinsic instability, as proposed in theoretical work^^ ^l^. 
Moreover, our model implies a continuous variation of 
the charge on the Fe ions (between 2+ and 3+, although 
the amplitude is probably smaller) clearly demonstrat- 
ing the presence of an incommensurate CDW at 150 K 
- to our knowledge the first example of such a structure 
in a half-doped oxide. Individual layers do not carry 
any net charge, and the local dipole moments cancel out 
over the whole bilayer, indicating that neither the ferro- 
electric nor the stacked antiferroelectric models can be 
correct. Rather, YbFe2 04 is an incommensurately mod- 
ulated antiferroelectric, as originally proposed by Ya- 
mada. We suggest, therefore, that the electrical data 
should be critically re-examined. It is well known^^ that 
the activated de-pinning of CDWs leads to pronounced 
non-Ohmic behaviour (observed in various organic and 
inorganic systemP^t^. The recent transport measure- 
ments performed on LuFe204'^ show precisely such an 
effect which, given the absence of breakdown of CO in 
an applied electric field^^, provides strong evidence for 
the existence of such a CDW state in this system. The 



observed pyroelectric currents^ could also be an arte- 
fact, since metastable states are notoriously ubiquitous 
in electric-field-driven CDWs^^. 

The phenomenology of YbFe204 above T* is equally 
intriguing: we have shown that twisting helices are still 
present and maintain a high degree of transverse coher- 
ence in the ab plane - at least 70 A (20a) - at all T. This 
diffuse scattering was previously interpreted in terms of 
two-dimensional orderin^^ however, the observed dif- 
fuse scattering is inconsistent with either intra-layer or 
intra-bilayer two-dimensional ordering, since these would 
produce rods, not helices of diffuse scattering. Significant 
c-axis correlations, much greater than the interbilayer 
distance, must be present in order to produce helices. 
As previously observed by Yamada, we have shown that 
the isointensity surfaces appear to follow the contour of 
the iso-energetic surfaces for modes of given wavevectors, 
as predicted by the random-phase approximation theory 
for critical fluctuations^^. In real space, this would result 
from a progressive population of modes within the quasi- 
degenerate helical manifold, starting from those that are 
closest to the ground state, so that the local coherence is 
maintained whilst allowing the system to explore states 
around the minima. We note that this scenario is very 
different from conventional CO transitions, which are 
driven by fluctuations of the amplitude of a single or- 
der parameter. 
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